Euclidean random matrix theory: low-frequency non-analyticities 

and Rayleigh scattering 

Carl Ganter 

Institut fiir Rontgendiagnostik, Technische Universitdt Miinchen, 

o ; 

, Ismaninger Str. 22, D-81675 Miinchen 

■ Walter Scliirmacher 

^ Institut fiir Physik, Universitdt Mainz, D-55099 Mainz, Germany and 

Physik- Department E13, Technische Universitdt Miinchen, D-85747, Garching, Germany 

Abstract 

By calculating all terms of the high-density expansion of the euclidean random matrix theory (up to 
' second-order in the inverse density) for the vibrational spectrum of a topologically disordered system 
we show that the low-frequency behavior of the self energy is given by Ti{k,z) oc iP'z'^l'^ and not 

claimed previously. This implies the presence of Rayleigh scattering and 
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long-time tails of the velocity autocorrelation function of the analogous diffusion problem of the form 



PACS numbers: 65.60. +a 



I. INTRODUCTION 

Rayleigh scattering i.e. the fact that the mean-free path of weakly scattered waves varies 
as cu^^'^^^^ in a ci- dimensional disordered medium as w — )■ 0, is widely believed to be a general 

n 

property of quenched disordered matter. However, recently it has been claimed [3[ that a 
harmonic system with displacements Ui{t) obeying 

d^ 

= - ^tij{ui{t) - Uj{t)) , (1) 

j 

where i and j denote random sites rj j in d-dimensional space, would have wave-like excitations, 
which have a line- width (inverse mean-free path), varying with u"^ instead of u'^ in d = 3. tij are 
force constants, divided the mass at the node i, which are assumed to depend on the distance, 
i.e. tij = t{rij). The claim of absence of Rayleigh scattering had been substantiated by a 
high-density expansion and a diagrammatical analysis This claim is not only astonishing 
with respect to the mentioned general view on waves in disordered media, but it is also in 
contradiction with the known analytic properties of the analogous diffusion system. If one 
replaces the double time derivative in ([T]) by a single one, one obtains the equation of a random 
walk among the sites 

^n,it) = ~J2^,,{n,it)-n,it)), (2) 
j 

where ni{t) give the odds for the walker to be at i at time t and = t{rij) is the hopping 
probability per unit time. Eq. ([2]) describes e.g. the motion of electrons hopping among 
shallow impurities in a semiconductor 5, 7, 8|. Such a random walk is known [6] to exhibit 
a long-time tail of the velocity- autocorrelation function (VAF) varying as Z(t) oc for 
t — > oo {gI, a feature shared with Lorentz models jsH]. In fact, the Laplace transform 
of the VAF is the frequ ence -dependent diffusivity D{z = iu + e), which has, according to 
the Tauberian theorems [,12j a low- frequency singularity D(z) — )■ z'^^'^, \z\ — )■ 0. Now, in the 
analogous vibrational problem this quantity corresponds to the square of a frequency- dependent 
sound velocity D{z = —u"^ + ie) = v'^{z). The imaginary part v"{uj) of the latter is related to 
the mean-free part via 

i{u) \v{z)\^ ^ ' 

This gives i oc u~^^^^\ i.e. Rayleigh scattering. We conclude that the long-time tail of the 
VAF in the diffusion problem is mathematically equivalent to the Rayleigh-scattering property. 
In the following we set the dimensionality d = 3 and calculate explicitly all irreducible diagrams 
(self-energy diagrams) up to second order in the inverse density p^^ = V/N, where is the 



number of sites and V the volume. We show that to this order the self energy is proportional to 
fc^z^/^, z = ico + e (diffusion) or z = — + ie (sound), and not as claimed in Refs. |i3|] oc k'^z^^'^. 
We also show, why the so-called cactus approximation for a self-consistent theory erroneously 
leads to a non-analyticity z^^"^ instead of z^^"^. 



II. FORMALISM 



As in refs. pi] we start from a high-frequency (z) and high-density (p = N/V) expansion of the 
averaged propagator 



CO ^ ^ 



p=l io...ip 

Here M is a matrix with off-diagonal elements Mjj = and diagonal elements Ma = — X^^^j ^i^- 
t{k) = t(k) is the (i- dimensional Fourier transform of t{r). 

The configurationally averaged Green's function can now be expressed in terms of the irreducible 
self energy T, {k,z) as follows 

^^'^'^^ ^ z- p[t{k) -t{0)]-J: {k, z) ^ z + D{z)k^ 
The frequency-dependent diffusivity /sound velocity is given by 

1 92 



D(z) =v(z? 



t{k) + i:{k,z) 



(5) 

fc-5-O 



2dk 

For simplicity, we assume complete site disorder (i.e. the radial pair correlation function 
g (r) = 1). Therefore t (k) is simply the Fourier transform of the transition rate t (r) jl3| . 
We denote the unrenormalized part of the Green's function by Go- 

'■= z-p[til)-tiO)] 



In analogy to the approach in Ref. [3|, it is helpful to consider a high-density expansion of the 
propagator, which is in turn determined by an analogous expansion the self energy S {k, z) 

00 

S(A;,^)=:5^p-"S(")(A;,;2), (7) 

n 

As outlined in [3|] , the index n counts repetitions of sites in the high frequency / high density 



In the following, we will derive exact results for n = 1 and n — 2. 

To this end, we will use diagrammatic representations, to distinguish topologically different 
contributions to the self energy. 

• Off-diagonal matrix elements, associated with a site change. Open circles will always 
indicate start and end point of a diagram. 

t1, = t (ri2) 1 O— 2, (8) 

• Part of the diagonal element, associated with site 2 - its value is determined by the sum 
rule. Since there is no site change, only one open circle is drawn. 

4 = -t(ri2):= 1 0- ► 2 (9) 

• The unrenormalized propagator 

Note that the propagator contains a diagonal part (sites 1 and 2 are the same, i.e. the diagram 
has length zero) 

formally obtained as Qo (z) = lim^^oo Co {k, z). In most cases, this requires no special attention. 
Exceptions, when these terms need to be explicitly excluded to preserve irreducibility will be 
mentioned below. 

III. FIRST-ORDER DIAGRAMS: S^^) (it, z) 

This case is comparably trivial and requires the addition of four diagrams, since the first and 
last connection can refer to an off-diagonal (O) or a diagonal (D) transition rate. 

00: 




OD: 




dp 

(2^ 



t (k - p) t (p) Go (p, z) 



DO: 




dp 

(2^ 



t (k - p) t (p) Go (p, z) 



DD: 



/ \ 
/ \ 

/ V 



p I -^t'{k-p) Go{p,z) 



Added together: 



p I -^[t{k-p)-t{p)fGo{p,z) 
(27r) 

Since we are mainly interested in the imaginary part of 



lim lim S (k, z) 

Z-5-0 fc-5>0 

we will have to examine the bracket in Eq. f|T6l) in the limit — )■ 



l,m<(k-p)-t(p) = ^:M.kp + i'"(f)-f-''(P' 



We therefore obtain 



P 



[kpf + \^^-k' + 
2 p 



lim lim [t(k- p) = c - [kp]VC p^) 



with some usually nonzero constant c. With the additional factor from the threedimensional 
integral, we obtain from the diffusion pole of Gq z) 



IV. SECOND-ORDER DIAGRAMS: (fc, z) 

This case is considerably more complex. It turns out to be advantegeous to consider topologi- 
cally different groups of irreducible diagrams separately: 







(l*2---2---l 


V(2) 


{k,z): 


(12---2---1) 


ym 


{k,z): 


(l*2---l---2 




{k,z): 


(12---1---2) 




{k,z): 





Unlike in ■ ■ ■ at least one additional site index (differnt from 1 and 2) needs to be contained in 



The complete and exact second-order self energy is then just the sum of these partial contribu- 
tions: 



S(2) {k, z) = S(2) {k, z) + (A;, z) + Sf {k, z) + (fc, z) + (fc, z) (21) 



A. (A;, z): Irreducible Diagrams (1 * 2 • • • 2 • • • 1) 

Here, 16 diagrams need to be distinguished jl^ : 

0000; 



2-i-O fe-i-O 



lim lim Im [ S^^^ ik, z) 1 oc z^^^k^ 



(20) 




(22) 



OOOD: 




-p'Jdpdqt{k- p) t' {p) t {q) Gl {p) Go (q) (23) 



OODO: 




f = -p" j dp e ip) t (p - q) Gl (p) Go (g) (24) 



OODD: 




/ Jdpdqt{k- p) ^2 (p) i (P - q) (P) Go (g) (25) 



ODOO: 




j dp ciq (p) t (p - q) t (g) (p) Go (g) (26) 



ODOD: 



j \ = j d^pdcitili-^) tip) t{p-c)t{q) Clip) Go{q) (27) 



ODDO: 




ODDD: 




= j dq e (p) e (p - q) Gl ip) Go (g) (28) 



p'J dpdqtik - p) t{p) t' (p - a) Clip) Go{q) (29) 




jdpdqt(k- p) (p) i (q) Gl (p) Go (q) (30) 



DOOD: 




DODO: 




DODD: 




j dv dci e (k - p) t (p) t (q) Gl (p) Go (q) (31) 



) = p'Jdpdcit{k-p)t'{p)t{p-ci)Glip)Goiq) (32) 



p' J dp dq e (k - p) t (p) t (p - q) (p) Go (g) (33) 



DDOO: 



i 

I 
I 



Jdpdqt{k-p)t (p) t (p - q) t (g) (p) Go (q) (34) 



DDOD: 





Si^) (fc, z) = (38) 
p'J dpdci[t{k-p)-t{p)f [t(p-q)-t(p)] [t(p_q)-t(g)] dip) Go{q) 

Two cases need to be distinguished: 

p is small: Because of Eq. fll9l) . the first squared bracket delivers a factor p^. Additional 
factors result from the third bracket and the 3D integration, respectively, so that 
we finally arrive at a p^ factor. From the identity {p, z) oc -^Gq (p, z) , we obtain a 
nonanalyticity oc z^^"^. Note that uneven occurrences of p and/or q, such as an isolated 



q is small: Now the second bracket delivers an additional nonanalyticity [15] g^, which again 
produces a z^/"^ nonanalyticity. 



We therefore conclude for this group of diagrams 

lim lim Im [ S^^) (A;, z) 1 oc z^''^k'^ (39) 



B. S^^^ {k, z): Irreducible Diagrams (1 2 • • • 2 • • • 1) 



Since sites 1 and 2 are connected directly via t'j'g or t^2 5 the diagrams contain only three t's and 
two Go's. Therefore, only eight diagrams are possible: 



000; 





ODD: 



t I. 



p dpdqt{k-p)t ip) t (p - q) Go (p) Go (g) (43) 




pjdpdqt{p)t{p-q)t (g) Go (p) Go (g) (44) 



DOD: 




^ p j c^pdqt(k-p) i(p-q) i(?) Go(p) Go (g) (45) 



DDO: 








p / dp dq t (p) e (p - q) Go (p) Go («) (46) 



DDD: 



-p / (ip<iq«(k-p) «2(p-q) Go(p) Go(?) (47) 



Added together: 



P 



^fik,z)= (48) 

i y dpdq[t(k-p)-t(p)] [t(p-q)-t(p)] [t{p-q)-t{q)]Go{p) Go{q) 

Because of f|T8l) . the first bracket gives us only a /c^. Depending on whether p or g are small, 
the third, respectively second, bracket delivers the required additional p^, respectively g^, to 
obtain: 



lim lim Im 

2-5>0 fc-s-0 



(A;,z)l ocz3/2fc2 (49) 



Note a particular property of diagrams DDO and DDD: 

If the propagator in the middle of the diagram collapses to the diagonal Qo, as explained above, 
these diagrams are f actor izable at site 1 and therefore not irreducible anymore. To avoid double 
counting, the diagonal term Qo must therefore be subtracted from this propagator. 
It can be easily verified, however, that the 2;^/^-nonanalyticity is not affected by this. 

(2) 

C. {k, z): Irreducible Diagrams (1 * 2 • • • 1 • • • 2) 

The crossover topology of these 16 diagrams leads to more intricate convolution integrals. 
Based on the following abbreviations 

a := k — p — q 

(3 ■■= P 
7 := q 

(50) 

a := p + q 



d(3 d-y Go (a) Go (/3) Go (7) ■ t{cx)t {(3) (7) 




dl3 dj Go (a) Go (13) Go (7) • (a) t (a) t (7) 



OODD: 



ODOO: 



ODOD: 




(54) 



(55) 



(56) 



df3 d-f Go (a) Go (/3) Go (7) ■ (a) t (/3) t (7) 



ODDO: 



ODDD: 



DOOO: 




(57) 



(58) 



(59) 



- p2 / dl3d-i Go (a) Go (/3) Go (7) ■ t{a)t (a) t (/3) t (7) 




J j d(3 dj Go (a) Go (/3) Go (7) ■ t{cx)t (b) t (c) t {(3) 




d(3 dj Go (a) Go (/3) Go (7) • t (a) (a) t (7) 




P 



J Jd/3dj Go (a) Go Go (7) ■ (b) t (c) t (/3) 




dl3 dj Go (a) Go (/3) Go (7) ■ t{cx)t (b) t (c) t (7) 





dp dj Go (a) Go Go (7) ■ (b) t (c) t (7) 



DDDD: 




(66) 



Collecting Terms 

To derive a usable expression for the sum of these 16 diagrams, we have to exploit the symme- 
tries of the problem. 

Transforming variables p, q — ^ p, q allows to arbitrarily permute a, /3, 7 under the boundary 
condition that a, b, c perform the same permutation, as indicated in the following table: 

a /3 7 a b c 

a 7 /3 a c b 

/3 a 7 b a c 

/3 7 a b c a 

7 a /3 c a b 

7 /3 a c b a 

For example, the transformation 



p := k — p — q und q P (68) 

leads to 



(67) 



k - p - q ^ p (69) 

p ^ q (70) 

q ^ k-p-q (71) 

p + q k - p (72) 

k-p ^ k-q (73) 

k - q p + q (74) 

and thus the permutation a (3^ ahc — )■ (3^ ahca. 

The product Go (ct) Go (/3) Go (7) is invariant with respect to these permutations and because 
of 




docdfi^ j j dad-f^ J j djSd-y (75) 

the integration variables can be chosen freely. After suitably regrouping the t-factors, we obtain 
for the sum of all diagrams the expression 

K := 0000+ +BBBT> 

^ J Jdf3dj Go (a) Go {/3) Go (7) x 

X [t{f3)-t{h)]-[t{j)-t{c)]-[t{j)-t{f3)]-[t{a)+tm 

Since we integrate over two variables only, we have to eliminate one Go-factor. To this end, we 
apply a partial fraction decomposition 

p [t (7) -tm-Go ((3) Go (7) = [Go (7) - Co (/3) ] (76) 
Some further permutations and regrouping lead to 

K = p2 y y ^^^^ Go (/3) Go (7) X 

X [t{a)-t{a)]-[t{cx)-t{j)]-[t{j)-t{c)] 
After reinserting the above definitions, we can expand for small k: 

i(a)-t(a) = t(k-p-q)-i(p + q) 



t(c)-t(7) = t(k-q)-t(g) 
q 

We thus obtain linijt^o-f^ oc /c^, as required. Since we only consider the lowest-order term, we 
may use 



i(a)-i(7)=i(a)-t(7) + C»(fc) 



(77) 



and set 



t{ot)-t{^)^t{p + <l)-t{q) 



(78) 



We thus finally arrive at 



K = J J dpdq G'o(p)Go(q) x 
t' (p + q 



(79) 



X 



|p + q 



■ k (p + q) 



t (p + q) - 1 (g) ; 



kq 



+0 {k^) 



A 



C 



For the analytical properties, we again have to consider two cases: 



q is small: 



Bracket C gives us a factor q, but brackets A and B approach a finite value for g — )• 0. 
But: After setting q = in A and B, both the integral over p (because of the kp term in A) 
and over q (due to the kq in C) vanish due to lack of rotational invariance. Consequently, 
we have to expand the fraction in A to first order in q, which provides us with an additional 
factor pq oc q. Now rotational invariance is preserved and a z^^"^ nonanalyticity is obtained. 



p is small: 

Bracket C remains finite and B leads to a factor pq oc p. 
Setting p = in ^4 leads to the following structure 



dpdq /(g) • [pq] • [kq]^ 



(80) 



In order to restore rotational invariance, we again have to expand in A to first order in p, which 
yields an additional pq oc p and leads to a z^/^ nonanalyticity. 
We thus conclude that 



hm hm Im [ S^^) {k, z) 1 oc z^^^k"^ 



holds. 



D. T,f^ {k, z): Irreducible Diagrams (1 2 • • • 1 ■ • • 2) 



Similar to above, 8 diagrams need to be considered: 



000; 



(81) 




= p J dpdqt{k-p + q)t{p)t (q) Go (p) Go (q) (82) 



OOD: 




,L ^> = -p [dpdqt{k-p + q)t{p)t{p-q)Go{p)Go{q) m 



ODD: 




DOD: 





DDD: 



t 



\ 



dp dq e (k - p) t (p - q) Go (p) Go {q) (89) 



Added together: 



(A;, z) oc (90) 
p J ciprfq[t(k-p)-t(p)][t(k-p + q)-t(k-p)] [t(p-q)-t(g)] Go{p) Go (g) 

The first bracket is proportional to k^. Depending on whether p or g is small, the third or 
second bracket yields the required additional or g^, respectively. 
Thus we have here too: 



lim lim Im 



{k,z) 



oc z 



3/2^2 



(91) 



E. si^^ {k, z): Irreducible Diagrams (1 ■ ■ • 1 ■ ■ ■ 1) 



Here only 4 cases are possible. 



00; 



-p Jdpdqt (p) t (p - q) t (q) Go (p) Go (q) 



(92) 



OD: 



DO: 




DD: 




-p y"c/pdq t (k - p) t (k - q) t (p - q) Go (p) Go (q) (95) 



Added together: 



{k, z) oc (96) 
p y dpdq [t(k-p)-t(p)][t(k-q)-t(g)] t (p - q) Go (p) Go (g) 

For order O {k'^), because of (|T8l) . the first two brackets muts be proportional to oc kp and kq, 
respectively. For rotational invariance, the last bracket needs to be expanded and thus provides 
the additional required pq term, to obtain: 



lim lim Im [ S^^) (k, z) 1 oc z^l'^k^ 

2^0 fc^O Lev /J 



(97) 



With fl2T|) it follows immediately 

lim lim Im [ S^^^ (k, z) 1 oc z^l'^k^ (9^ 



2-i>0 



V. CONCLUSION 

Working out term by term in the second-order self energy we have convinced ourselves that to 
this order the nonalytic behavior (1981) . which both leads to Rayleigh-type sound attenuation and 
to the correct long-time tail in the analogous diffusion problem is recovered. We have done this 
for = 3, but the generalization to any > 1 is straightforward. Our result is in contradiction 
to the claims in the publications [3]. In these publications a self-consistent equation for the 
self-energy is advocated, which consists in making the first-order result f|T6|) self-consistent, i.e. 
replacing the 0-th-order Green's function by the full one. Now, in performing a high-density 
expansion of this equation one easily convinces oneself that the corresponding diagrams are 

• the entire sum (A;, z); 

• the diagrams 000, OOD, ODO, ODD of ^ '{k,z), but not the remaining four dia- 
grams; 

• the diagrams DOO, DOD, DDO, DDD of J:f\k,z), but not the remaining four dia- 
grams. 

As the partial sums do not give the correct analytic properties, this is the reason, why the 



self-consistent scheme advocated by js] does not lead to Rayleigh scattering. We shall publish 
shortly a self-consistent scheme, which includes Rayleigh scattering. 
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